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Abstract. The McLerran-Venugopalan model provides a framework which allows one 
to compute the gluon distribution function of a very large nucleus from the equations of 
QCD, provided that the longitudinal momentum fraction, x^, is sufRciently small. The 
source of color charge for this computation may be thought of as a crepe moving along 
the z axis at the speed of light. We refine the MV model by allowing for the presence 
\^ ' of non-trivial longitudinal correlations between the color charges that comprise the 

. nucleons. We find that a consistent treatment forces us to consider a pancake-like 

source which moves at slightly less than the speed of light. Our calculation allows us 
. to consider larger values of than were allowed in the original MV model. 
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Several years ago, McLerran and Venugopalan realized that for large enough 
CL|| nuclei at small enough values of the longitudinal momentum fraction Xp , it ought 

Q^' to be possible to compute the gluon distribution function using QCD [1]. Based 

1^ ■ on this observation, the framework known as the McLerran-Venugopalan model 

(MV model) was subsequently developed [2-5]. Recently, it was shown [6] that 
the infrared divergences present in the MV model may be cured by capturing the 
^ ■ physics of confinement via a color-neutrality condition to be imposed on the charge 

density correlation function used as input to the MV model. In this talk, I will 
describe work [7] on extending the MV model to larger values of x^. 

The MV model as originally formulated in Refs. [1-5], is restricted to very small 
longitudinal momentum fractions Xp ^ A^^/^/(ma), where A is the number of nu- 
cleons, m is the nucleon mass, and a ~ AqJ^j-, is the nucleon radius. In this regime, 
the longitudinal resolution of the gluons is so poor that they probe distances which 
are much longer than the (Lorentz-contracted) thickness of the nucleus. Thus, all 
of the quarks inside the nucleus effectively have the "same" value of the longitu- 
dinal coordinate x~ . At each value of the transverse coordinate x describing this 
crepe-shaped {i.e. very thin) nucleus, the color charges from a large number of 
valence quarks must be summed, resulting in a color charge density which is in a 
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high-dimensional representation of the gauge group. This is a necessary condition 
for a classical treatment to be vahd. In going to larger values of x^, we find that the 
longitudinal resolution of the gluons improves, and we begin to see the longitudinal 
structure of the nucleus. In order to include this longitudinal structure, we are nat- 
urally led to a pancake-shaped geometry [i.e. one with a finite non-zero thickness). 
The details of this fully 3-dimensional calculation are contained in Ref . [7] . 

For a classical treatment to be valid, not only must the color charge be in a 
large representation of the gauge group, but the gauge coupling should also be 
weak. McLerran and Venugopalan [1] argue that the running coupling ought to be 
evaluated at the scale /x^, which is set by the charge- squared per unit transverse 
area. For large enough nuclei, //^ » -^qcd) implying that Q;s(/i^) <C 1. 

Thus, we begin our computation of the gluon distribution function for a large 
nucleus by solving the classical Yang-Mills equations describing a pancake-shaped 
distribution of color charge moving along the z axis at nearly the speed of light. 
The result is a non-linear expression for the vector potential A{x~,x) in terms 
of the charge density p{x~,x). In the spirit of the Weizsacker- Williams approx- 
imation [8], we extract the gluon number density from the two-point correlation 
function, {A{x~ , x)A{x'~ , x')) . We replace the quantum mechanical average im- 
plied by the angled brackets with a classical average over an ensemble of nuclei. 
This ensemble is specified by inputting the two-point charge-density correlator 
(pp) ~ V{x^,x). Furthermore, we assume that the correlations are Gaussian. 
Confinement is incorporated into the calculation at this stage in the form of a color 
neutrahty condition on V [6]. When V satisfies the color neutrality condition, the 
two-point correlation function is infrared finite, and may be Fourier-transformed 
to momentum space, producing a gluon number density dN/dXpd?q which is dif- 
ferential not only in x^, but in the transverse momentum q as well. 

In the limit A"^/^ <C 1, the result of this rather lengthy calculation reads^ 

Although complicated in appearance, Eq. (1) is made up of several easily- 
understood parts. The prefactor shows that at lowest order the number of gluons 
is simply proportional to the number of quarks. The lowest order result is governed 
by the function 

^l^., - 2 y 4^2 e [g2 + (^^^)2]2' 

and is what would be obtained by considering an Abelian theory. The non-Abelian 
corrections to this result are contained in the exponential factor, which depends on 
two quantities. First, the the spatial dependence is determined by 



^) Eq. (1) has been written assuming cylindrical geometry for the nucleus. For a full discussion 
of the (rather weak) geometric dependence of the result, see Ref. [7]. 
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The strength of the non-Abehan corrections is set by the factor 

= ik^W^ ~ 8™;^'"A?,co. (4) 

These effects are most prominent in very large nuclei. For uranium we have A'oo ~ 5 
or BAqqq. To obtain X^^ — 20Aq(-.j5 (as is employed in the plots below) requires 
of order 10^ nucleons. Finally, we note that if we set = in Eqs. (l)-(3), we 
explicitly reproduce the original MV result [5]. 

The functions appearing in Eqs. (2) and (3) are not finite at qr = unless the 
charge density correlator V satisfies the requirement of color neutrality. As ex- 
plained in Ref. [6], a key consequence of confinement is the appearance of color 
neutral nucleons. Mathematically, this consequence may be implemented as a con- 
straint on T>: 

J dA-(fA V(A-, A) = 0, or V(0, 0) = 0. (5) 

Eq. (5) ensures that widely-separated nucleons are essentially uncorrelated. In 
order for (5) to be true, V must contain a length scale. Not surprisingly, this scale 
turns out to be the nucleon radius a ~ Aq^j-,. When combined with the assumption 
of rotational symmetry in the transverse plane, Eq. (5) is sufficient to render the 
functions contributing to the gluon number density completely infrared finite: the 
integrals get cut off at g ~ ~ Aqcd- 

Aside from the neutrahty condition (5), the correlation function T> is unspecified. 
In order to illustrate the general features of the gluon number density (1), it is 
convenient to choose the following form for T>: 

'^{^f,q) = ^ - 7 — so ^ / ^- (6) 

In the context of Kovchegov's nuclear model [9], the correlation function given in 
Eq. (6) corresponds to quarks which are distributed within the nucleons according 
to the weight 



At large distances, the modified Bessel function produces a Yukawa-like behavior 
in this probability distribution. 

Fig. 1 exhibits the behavior of the fully differential gluon number density as a 

function of the transverse momentum . We see that there is saturation: the 
number of gluons increases as ^ to some maximum value and then stops 



growing. Qualitatively, these distributions are very much like those obtained by 
MuUer from the point of view of onium scattering [10,11]. 

The gluon distribution function resolved at the scale is related to the fully 
differential gluon number density by 
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In Fig. 2 we multiply the fully differential distribution by q^: on the semi-log scale 
used this produces a "true" representation of where the important contributions 
to Eq. (8) are. We see that the very low region does not play a significant 
role provided is not too small. From the plot we see that for oo, the 

gluon distribution function becomes insensitive to the presence or absence of the 
non- Abelian corrections. In fact, we can prove under rather general circumstances 
that 
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independent of V [6]. 

In Fig. 3, we exhibit XpgA^x^, 
of Q^. At low values of Q^, the addition of the non- Abelian corrections reduces 
the number of gluons from the Abelian result, while by the time = 2500Aq(J£, is 
reached, the effect of the non- Abelian terms is negligible. Although we have plotted 



as a function of a;„ for several different values 
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FIGURE 1. Fully differential gluon number 
density plotted versus at = 0.0 and 0.1. 
The three curves are for X^o — (Abelian), 
5Aq(-.j-,, and 20Aqq]3. The gluon density satu- 
rates at small q^. 
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FIGURE 2. Fully differential gluon number 
density at =0.1 multiplied by q^ versus 
for Xoo = (Abehan), SAqj;;^, and 20Aq(.j3. 
The area under each curve accurately reflects 
the contributions to gA(xp ,Q^)- 



out to Xp — 1, the calculation is not reliable beyond ~ 0.25 [7]. For small x^ 
the gluon distribution function exhibits the pure 1/Xj, dependence characteristic of 
the original MV model. 

Finally, we present Fig. 4, which illustrates the deviation of the gluon distribution 
function from the naive expectation that for A nucleons we should obtain A times 
the result for a single nucleon. We see that at low values of the distribution 
function grows less rapidly than the number of nucleons as ^4 is increased, whereas 
at large Q^, the simple scaling expectation holds. 
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FIGURE 3. The gluon structure function 
XpQ^ {Xp , Q^) plotted versus x^. for several val- 
ues of Q^. The dashed curves illustrate the 

effect of ignoring the non- Abelian corrections, 
which are included in the solid curves. 
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FIGURE 4. Nuclear dependence of the gluon 
structure function Xpg^{Xp,Q^). For suffi- 
ciently small the number of gluons grows 
more slowly than the number of nucleons as A 
is increased. 
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